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Pulse  Combustion  Rockets  for  Space  Propulsion  Applications* 1 

Edward  B.  Coy* 

Air  Force  Research  Laboratory,  Edwards  AFB ,  CA  93524-7660 


Abstract 

Pulse  combustion  propulsion  devices 
are  currently  being  considered  as  alternatives  to 
conventional  constant-pressure  engines. 

Potential  advantages  include  reduction  or 
elimination  of  pumps  and/or  compressors,  and 
improved  Isp  for  a  given  feed  system  supply 
pressure.  In  this  paper  a  model  is  presented  for  a 
monopropellant-fueled,  constant- volume,  pulse 
combustor  which  includes  finite-rate  processes. 

A  zero-dimensional  model  is  used  for  the 
combustion  chamber  and  a  one-dimensional, 
quasi-steady  approximation  is  used  for  the 
nozzle  flow.  The  liquid  spray  is  assumed  to  have 
a  log-normal  distribution  of  spherical  droplets 
and  the  reaction  rate  is  based  on  a  strand  burner 
correlation.  This  model  was  developed  as  a  tool 
for  designing  an  experimental  rocket.  In  this 
paper  it  is  used  to  explore  the  time  and 
dimensional  scales  of  the  problem  and  to  predict 
the  performance  and  optimal  geometry.  The 
pulsed  propulsion  device  is  found  to  have  nearly 
identical  specific  impulse  as  the  steady-state  • 
engine  operating  with  the  same  mass  flow  and 
throat  area,  furthermore,  the  nozzle  optimizes  at 
the  same  area  ratio.  Pulsed  combustor  behavior 
is  found  to  depend  on  two  time  scales:  the  ratio 
of  the  heat  release  time  to  the  chamber  blow¬ 
down  time,  and  the  ratio  of  the  blow-down  time 
to  the  injector  pulsing  period.  Finally,  the  model 
is  used  to  assess  potential  benefits  of  pulsed 
engines  for  satellite  applications.  We  briefly 
consider  the  application  of  pulse  combustion 
devices  in  pressure-fed  satellite  propulsion 
systems  and  examine  the  effect  on  satellite 
mission. 

Introduction 

In  a  previous  paper  we  explored  the 
constant- volume  limit  of  pulsed  propulsion, 
where  the  combustion  chamber  was 
approximated  as  being  time  varying  but  spatially 
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uniform,  and  the  nozzle  flow  was  approximated 
as  one-dimensional  but  quasi-steady  .  The  heat 
release  process  was  assumed  to  completely 
consume  the  propellant  and  to  occur  without 
change  in  mass  in  the  combustion  chamber, 
either  because  the  injection  and  reaction  rates  are 
much  faster  than  the  rate  of  mass  flow  out  of  the 
nozzle,  or  because  a  hot-gas  valve  is  located  in 
the  exit  stream.  In  this  paper  we  extend  the 
analysis  by  considering  finite-rate  processes  such 
as  heat  release,  injection,  and  blowdown,  and 
inefficiencies  associated  with  incomplete 
combustion. 

The  type  of  pulse  combustor  considered 
here  operates  on  the  following  cycle.  Reactants 
are  injected  in  a  short  burst  during  the  low- 
pressure  portion  of  the  cycle.  The  spray  ignites 
when  exposed  to  hot  residual  products  of  the 
previous  cycle  and  the  chamber  pressure  rises  as 
the  decomposition  products  are  produced.  The 
chamber  pressure  reaches  a  peak  as  the  liquid  is 
completely  consumed  and  then  falls  as  the  gas 
exits  through  a  Delaval  nozzle,  generating  thrust. 

This  engine  cycle  may  have  advantages 
relative  to  the  pulsed  detonation  engine  cycle2  in 
that  it  may  be  possible  to  eliminate  the  ignition 
system  for  all  but  the  initial  starting  transient.  It 
may  also  be  possible  to  eliminate  the  purge  gas 
separating  products  and  reactants.  Finally,  it 
may  be  possible  to  fill  the  chamber  to  high  initial 
densities  without  generating  pressures  and  shock 
waves  characteristic  of  high-explosives  (109  Pa). 

To  achieve  pulsed  operation  the 
injection  time  must  be  less  than  the  heat  release 
time  which  in  turn  must  be  less  than  the 
blowdown  time.  Furthermore,  to  produce  large 
pressure  oscillations,  the  blowdown  time  must 
not  be  long  relative  to  the  injector  pulsing  time. 

If  this  order  is  not  maintained  then  the 
advantages  of  low-pressure  injection  cannot  be 
realized. 

Two  methods  of  injection  may  be 
considered:  fixed  frequency  and  fixed  pressure. 

In  general  fixed  frequency  is  preferred  because  it 
enables  precise  control  of  the  mass  flow  rate. 

Fixed  pressure  has  the  advantage  that  it  allows 
propellant  to  be  injected  at  the  point  in  the  cycle 
which  is  optimal  for  achieving  the  desired  heat 
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release  time.  Both  methods  result  in  stable 
operation,  therefore  from  an  analytical 
perspective,  there  is  no  difference.  In  this  paper, 
fixed  frequency  operation  is  assumed. 

Propulsion  devices  operating  on  the 
constant- volume,  or  “explosion-cycle”  have  been 
proposed  in  the  past.  Ref.  3  discusses  an  early 
attempt  at  building  an  air-breathing  device.  This 
device  was  fueled  by  a  gasoline  spray.  The  heat 
release  time  was  measured  to  be  6-9  msec  which 
produced  a  pressure  rise  of  approximately  a 
factor  of  two.  It  was  predicted  that  a  pressure 
rise  time  of  3  msec  was  needed  to  generate 
pressure  ratios  for  viable  operation. 

Detonation  wave  engines  can  complete 
reactions  on  the  order  of  1  msec  or  less  and 
generate  pressure  ratios  on  the  order  of  10  for 
air-breathing  engines.  For  rocket  engines 
operating  on  condensed  phase  propellants,  the 
pressure  ratios  could  potentially  reach  extremely 
high  values.  In  fact,  the  peak  pressures  must  be 
constrained  by  limiting  the  density  of  the  initial 
charge  to  avoid  spallation  and  extremely  heavy 
combustion  chambers.  Since  the  pressure  ratios 
must  be  limited,  it  is  worthwhile  considering 
alternative  means  for  achieving  combustion¬ 
generated  pressure  ratios.  In  this  paper  we  use  a 
correlation  developed  for  the  monopropellant 
nitromethane  to  show  that  very  large  pressure 
ratios  can  be  generated  from  the  decomposition 
of  a  spray. 


Model 

The  model  consists  of  standard 
equations  for  conservation  of  mass  and  energy  in 
a  control  volume.  A  “lumped-parameter” 
formulation  is  used.  The  uniform  chamber 
approximation  can  be  made  for  the  devices 
considered  in  this  study  because  the  transit' time 
for  a  pressure  wave  to  traverse  the  combustion 
chamber  is  much  less  than  the  time  scale  of 
pressure  changes  in  the  chamber.  A  typical 
value  for  speed  of  sound  is  1200  m/s  and  a 
typical  chamber  length  is  5  cm,  giving  a  transit 
time  of  4x1  O’5  seconds,  whereas  the  time  scale 
for  pressure  changes  due  to  combustion  and 
blowdown  are  on  the  order  of  10'3  seconds. 

The  control  volume  contains  the  gas  in 
the  combustion  chamber  up  to  nozzle  throat  but 
not  the  droplets.  For  the  continuous  processes  of 
droplet  burning  and  chamber  blowdown, 

mg  =  ms  -  me ,  (i) 


cLT  .  u*  .  n  ue 2 

ms  Cr~17  m‘  (h*  +^~me  (he  +  — )  (2) 
at  2  2 


+  Qcv  ~  Wcv 

where  the  subscript  g  denotes  gas,  s  denotes 
spray  and  e  denotes  exit.  The  velocity  of  flow 
from  the  droplets  is  assumed  to  be  negligible. 
The  energy  entering  the  control  volume  from  the 
spray  can  expressed  in- terms  of  an  adiabatic 
flame  temperature,  hs  =  CPTAFr .  The  exiting 
energy  can  be  expressed  in  terms  of  the  chamber 
stagnation  enthalpy,  cPTg  =  ht  +  u]  12  •  The  work 


term  accounts  for  expansion  of  the  control 
volume  as  droplets  in  the  chamber  are  consumed, 
Wcv  =  pg  ms !  ps .  Heat  transfer  is  neglected. 

Assuming  isentropic  flow  of  a  perfect 
gas,  the  exiting  mass  flow  is  calculated  from  the 
choked  flow  expression, 


m 


exit 


(3) 


where  T0  is  the  instantaneous  stagnation 
temperature  in  the  combustion  chamber.  The 
exiting  flow  is  assumed  to  be  choked  at  all  times. 

Mass  flow  from  the  droplets  is  obtained 
using  the  following  model.  The  initial  droplet 
size  distribution  for  each  injection  is  assumed  to 
be  adequately  described  by  the  log-normal 
distribution. 


n(D0)  = 


N 


ln2cr  D, 


xexp 


(ln£>0  -  In  Dg ) 
2  In2  cr 


(4) 


where  N  is  the  total  number  of  droplets,  og  is 


the  median  diameter  and  a  is  the  geometric 

standard  deviation.  As  the  droplets  are 
consumed  the  distribution  will  change.  In 
general,  burning  rate  is  dependent  on  diameter  as 
in  the  case  of  the  “D-squared  Law”  and  the 
distribution  transforms  as,  n(D)dD=n(D0)dD0, 

where  D  is  the  current  size  of  the  droplet  which 
was  initially  D0at  the  time  of  injection.  Thus  to 

obtain  n(D)  from  n(D0)  one  must  have 
dD0  jdD .  For  the  case  of  monopropellant 

burning  the  surface  regression  rate  is  the  same 
for  all  diameters  so  the  derivative  is  equal  to  1 
and  the  new  distribution  is  simply, 
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n(D,  D')  = 


-Jin  In1  a,  (. D  +  D ') 


(ln(g  +  P')-lnP£y 


jc  exp  - 


where  D'  is  the  amount  of  surface  regression 
which  has  occurred  since  the  time  of  injection. 
The  total  volume  of  liquid  remaining  is, 


V(D')  =  ]^—  n(D,  D")dD 

This  integral  has  the  following  solution, 


V(D')  =— | [f3(P')-3D' F2(D')+3D,2Fi(D')-D*  F0(D') 
6 


where, 

F„  (0')=4  «pf nlnDg  +  ^-\n2as 


x\\-erf\ 


\nD'-(lnDs  +nln2cr5  ) 
Voinov 


The  mass  flow  rate  from  liquid  to  gas  used  in 
eqns.  (1)  and  (2)  is  obtained  as  follows, 


dV  dU 

”h  ~PmmD  dD'  dt  • 


(9) 


The  first  derivative  on  the  RHS  is  obtained  from 
eqn.  (7).  The  second  is  obtained  from  an 
empirical  correlation  contained  in  ref(4). 


—[mis]  =-.632xl0"9  Phm  [Pa] .  (10) 

dt 


Equation  (7)  is  solved  numerically  along  with 
the  conservation  equations  for  mass  and  energy 
in  the  control  volume.  To  determine  the  liquid 
volume  and  its  rate  of  change  from  any  given 
injection,  the  current  value  of  D‘  is  subtracted 
from  the  value  at  time  of  injection. 

It  is  worth  noting  that  the  most  common 
method  for  modeling  droplet  burning  is  to 
discretize  the  distribution  into  several  size 
classes  and  represent  each  class  by  an  ODE.  If 
the  device  is  a  pulsating  combustor  then  each 
new  injection  adds  several  more  ODEs. 
Eventually  computing  speed  must  suffer.  The 
model  presented  here  requires  only  one  ODE  for 
any  number  of  injections.  Furthermore,  the 
results  are  an  exact  solution  of  the  assumed  log¬ 
normal  distribution. 


Combustion  inefficiency  occurs  when 
the  residence  time  of  the  propellant  is 
insufficient  for  complete  combustion  to  occur. 
Within  the  framework  of  the  0-D  or  lumped- 
parameter  model  for  the  combustion  chamber, 
the  residence  time  of  the  droplets  is  assumed  to 
equal  the  residence  time  of  the  gas.  Therefore, 
for  each  injection  we  can  write  the  following 
expression  for  the  rate  of  change  of  the  number 
of  droplets  remaining  in  the  chamber, 


n  mi 

The  same  approach  used  in  calculating 
the  extent  of  vaporization  can  be  used  here  to 
eliminate  the  need  to  solve  an  equation  for  each 
injection.  The  formal  solution  for  eq.  (1 1)  is, 


The  first  integral  in  (12)  is  solved 
numerically  along  with  the  conservation 
equations.  The  value  at  the  time  of  injection  is 
stored  as  the  second  integral.  When  the  number 
of  droplets  remaining  in  the  chamber  is  less  than 
0.01%  of  the  number  injected,  the  injection  is 
eliminated. 

The  injection  process  occurs  as  discrete 
events.  A  finite  duration  injection  must  be 
simulated  by  breaking  the  continuous  event  into 
several  discrete  pulses  spaced  over  the  interval. 
Since  the  injection  occurs  instantaneously,  mass 
flows  do  not  occur  and  the  mass  of  the  gas  is 
unchanged.  Therefore, 

A Ecv  =  -A Wcv .  Assuming  isentropic 
compression  the  work  performed  on  the  gas  can 
be  calculated  as, 

AWC,  =  J PdVm^L{vtr  -V^). 

J  |  -  y 

v.  • 


.  03) 

where  V2  =  -  V5  .  The  step  change  in  gas 

temperature  from  the  compression  that  occurs 
during  an  injection  is,  A Tg  =  A Wcv  lmgCv  . 

The  exiting  flow  is  assumed  to  be 
choked  at  all  times.  If  the  nozzle  exit  pressure 
falls  below  the  level  at  which  a  shock  can  enter 
the  nozzle  a  flag  is  set.  The  data  presented  in 
this  paper  is  for  cases  where  a  shock  did  not 
enter  the  nozzle. 

A  monopropellant  simulant  with  the 
following  properties  was  used  for  all  the  cases 
presented. 
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Gas  Molecular  Weight 

20  gm/mol 

Specific  Heat  Ratio 

1.2 

Density  of  Liquid 

1000  Kg/m3 

Adiabatic  Flame  Temperature 

2500  K 

The  instantaneous  thrust  is  calculated 
using  the  ideal  thrust  equation, 


T  =  PA' 


2  r 
jr- 1 
(14) 


_2_r 

O'  +  l, 


1- 


\ 

PJ 


+  {pe-pa)At 


The  minimum  exit  pressure  at  which  a  shock  will 
enter  the  nozzle  is  obtained  from  the  expression, 


P_ e_ 
P 


2  Y 
y+l 


//  -  \Mr- 1) 

•  (15) 


The  exit  Mach  number  is  a  function  of  the  area 
ratio. 


The  performance  parameter  of  interest 
is  the  specific  impulse  or  Isp.  For  a  steady  • 
thruster,  Isp  is  defined  as, 

Isp  =  J—,  (16) 

gc™ 

where  gc  is  the  gravitational  constant.  For  an 
unsteady  device  Isp  is, 

Isp  =  -  P.(t)dt  >  (17) 

gc  J  rn(t)dt 

where  the  integration  is  over  the  duration  of 
operation,  or  over  a  single  pulse  if  the  device  has 
reached  limit-cycle  conditions,  or  over  a  sample 
of  pulses  if  there  is  significant  pulse-to-pulse 
variability.  The  mass  flow  rate  as  defined  here 
includes  the  amount  of  propellant  that  is  required 
to  pump  fuel  into  the  chamber.  For  simplicity, 
we  will  assume  a  pressure  fed  system  so  that  this 
mass  need  not  be  explicitly  accounted  for. 

In  addition  to  the  performance,  we  are 
also  interested  in  the  magnitude  of  pressure 
oscillations, 

P  ~  0° CMAX  ~Pc,Wn)/PsS  ’ 

where  Pss  is  the  steady  state,  or  high  frequency 
limit  of  the  chamber  pressure  for  a  given  throat 
area  and  mass  flow  rate. 

Performance 

In  this  section  some  example 
calculations  are  presented  which  show  the  effects 


of  pulsed  operation  on  performance.  Cases  1-4 
in  Table  1  show  the  effect  of  injection  frequency 
for  fixed  geometry  and  average  mass  flow  rate. 
Case  1  shows  that  injecting  at  5000  Hz  results  in 
nearly  constant  chamber  pressure  (P*=0)  and 
represents  a  baseline  for  comparison.  Cases  2-4 
show  that  as  frequency  is  decreased  from  5000 
Hz  to  100  Hz,  the  Isp  decreases  slightly  from 
233.30s  to  232.85,  or  0.2%,  while  the  minimum 
chamber  pressure  decreases  from  2.51  MPa  to 
0.28  MPa.  Thus  under  the  conditions  shown 
here  it  is  possible  to  achieve  nearly  the  same 
performance  as  a  conventional  rocket  with  an 
injection  system  that  discharges  at  one  tenth  the 
pressure. 

Chamber  pressures  for  Cases  2-4  are 
plotted  in  Figure  1.  Three  complete  cycles  are 
shown.  There  is  a  trend  towards  larger  chamber 
pressure  oscillations  at  longer  pulse  periods. 
Decreasing  the  frequency  at  fixed  engine 
geometry  and  mass  flow,  decreases  the  pressure 
at  the  time  of  injection. 

In  cases  5  and  6  the  area  ratio  is  varied 
above  and  below  the  baseline  case  to  determine 
if  the  optimum  value  differs  from  that  of  the 
equivalent  steady  state  engine.  A  baseline  area 
ratio  of  4.26  was  used  in  Case  4  which  is  the 
optimized  value  for  a  steady  rocket  operating 
with  the  same  average  mass  flow  and  throat  area. 
In  both  cases,  Isp  decreases  as  Ae/A*  moves 
away  from  4.26  indicating  that  the  optimum 
value  is  approximately  the  same  as  that  of  the 
steady  engine. 

In  cases  7  and  8  the  performance  in 
vacuum  is  compared.  An  area  ratio  of  100  was 
arbitrarily  selected.  Once  again  the  performance 
penalty  of  pulsating  operation- is  approximately 
0.2% 

An  alternative  basis  for  comparing  the 
performance  of  steady  and  pulsed  systems  is 
constant  injection  pressure.  Again  using  Case  4 
as  the  baseline,  in  Case  20  the  nozzle  throat  area 
was  decreased  while  average  mass  flow  rate  and 
injection  frequency  were  held  constant. 

Chamber  volume  was  also  decreased  to  maintain 
approximately  the  same  value  of  P*.  The  nozzle 
area  ratio  was  optimized  for  the  new  value  of 
Pss*  The  Isp  of  Case  20  is  273  s,  an  increase  of 
17%  relative  to  the  baseline  Case  4.  The 
increased  performance  is  due  to  the  higher  value 
of  PSs  ,■  15.3  MPa  for  Case  20  versus  2.51  MPa 
for  Case  4.  Case  20  assumes  an  ambient 
pressure  of  1  bar.  Figure  2  shows  additional 
comparisons  for  a  range  of  ambient  pressures 
from  1  bar  to  vacuum.  Nozzle  area  ratios  were 
optimized  based  on  matching  exit  pressure  to 
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Table  2  Pulse  combustor  performance  calculations 


ambient  assuming  a  chamber  pressure  of  PSs-  If 
the  required  area  ratio  exceeded  100  then  100 
was  used.  Clearly  the  benefit  to  performance  of 
increased  Pss  diminishes  as  ambient  pressure  is 
reduced.  This  is  simply  a  consequence  of  the 
diminishing  effect  of  increasing  pressure  ratio. 

The  above  discussion  shows  that  pulsed 
propulsion  has  essentially  the  same  performance 
as  steady  for  equal  values' of  Pss*  Furthermore, 
nozzle  area  ratios  optimize  at  the  same  value. 
The  potential  advantage  of  pulsed  propulsion  lies 
in  the  possibility  of  reducing  feed  system 
pressure  for  a  given  value  of  Pss-  This  may 
result  in  feed  system  components  which  are 
lighter,  less  expensive  and  more  robust. 


Pulsed  Propulsion  Scaling 

In  this  section  some  example 
calculations  are  given  to  illustrate  scaling 
principles. 

If  it  is  assumed  that  PSs  is  fixed  and  that 
the  heat  release  time  is  very  short,  then  the 
injection  period  and  chamber  volume  are  the 
only  independent  variables.  In  figure  1  the  effect 
of  increasing  the  injection  period  at  fixed  volume 
was  shown.  It  is  also  possible  to  show  that 
reducing  the  chamber  volume  at  fixed  injection 


period  has  an  identical  effect.  In  figure  3  the 
volume  and  injection  period  are  changed 
together.  The  baseline  Case  4  is  plotted  with 
Case  9,  in  which  pulse  period  and  chamber 
volume  are  doubled,  and  Case  10,  in  which  they 
are  halved.  In  each  case  the  chamber  pressure 
oscillates  between  the  same  limits.  The  injection 
pressure  and  performance  in  each  case  are  the 
same,  but  the  cases  with  short  injection  periods 
have  smaller  chamber  volumes. 

Changing  the  chamber  volume  affects 
the  characteristic  blow  down  time,  tBD  .  An 
estimate  for  tBD  time  can  be  defined  as 

tBD  =  rnc  ! mexit  =  V  I  A*  ^R  TC  f(y)  ■  If  we 

consider  the  ^ ~T dependence  as  a  weak 
function  of  the  chamber  pressure  then  for  the  gas 
considered  here  tBD  ~  1.5x10~3V7A\  Since  tBD 
is  a  linear  function  of  the  chamber  volume,  then 
the  effect  of  halving  or  doubling  volume  in 
Cases  9  and  10  was  to  halve  or  double  tBD  .  For 
the  conditions  given  for  figure  3,  we  can 
conclude  that  chamber  pressure  is  a  function  of 

*BD  ^ pulse  only- 

This  result  can  be  extended  to  include 
any  value  of  PSs  if  the  chamber  pressure  is  scaled 
by  Pss-  In  Case  1 1  the  throat  area  and  the 
injection  period  are  halved  while  in  Case  12  both 
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Table  3  Pulse  combustor  performance  -  continued 


are  doubled,  leaving  tBD  / tPULSE  unchanged. 

The  average  mass  flows  are  the  same  in  each 
case  so  the  changes  in  throat  area  result  in 
different  values  for  P ss.  If  time  is  also  scaled  by 
tPULSE  then  Figure  4  shows  that  the  non- 
dimensionalized  chamber  pressures  traces  are 
identical.  The  slight  scatter  observable  in  the 
traces  is  a  result  of  a  small  but  finite  heat  release 
times.  Figure  5  extends  the  results  of  figure  4  to 
show  the  effect  of  varying  tBD  I tPULSE  . 
Decreasing  tBDftPULSE  results  in  larger  relative 
changes  in  chamber  pressure. 

Figure  6  shows  the  effect  of  heat  release 
time,  tHR  .  There  is  no  simple  analytical  estimate 
for  tHR  so  the  definition  was  arbitrarily  chosen  to 

be  the  time  required  for  95%  of  the  injected  mass 
to  be  consumed.  The  most  direct  way  to  effect 
the  heat  release  time  is  through  the  droplet  size 
distribution.  Case  13  is  the  same  as  4,  but  Dg 

has  been  increased  to  10  micrometers.  This 
results  in  a  value  tHR  /tBD  =  0.55 .  The  extent  of 


chamber  pressure  oscillations  as  represented  by 
P*  is  reduced  from  3  to  1.7.  In  Case  14  Dg  =  10 

but  crg  is  increased  to  2.0  so  the  distribution  is 

no  longer  monodisperse.  This  results  in  further 
attenuation  of  chamber  pressure  oscillations. 

The  larger,  slower  burning  droplets  introduce  a 
time  delay  in  the  heat  release  which  dampens 
chamber  pressure  oscillations.  The  magnitude  of 
chamber  pressure  oscillations  can  be  recovered 
by  reducing  tBDftPULSE . 

Rapid  heat  release  is  desirable  because 
it  allows  for  a  smaller  value  for  tBD  which 

translates  into  a  smaller,  lighter  combustion 
chamber.  However,  if  heat  release  is  too  rapid 
then  it  will  not  be  possible  to  complete  injection 
during  the  low-pressure  part  of  the  cycle.  Figure 
7  shows  the  effect  of  finite  duration  injection  on 
the  rate  of  chamber  pressure  rise.  The  injection 
event  is  broken  up  into  10  discrete  pulses  evenly 
spaced  in  time.  Case  4  is  used  as  the  baseline, 
but  droplet  size  has  been  increased  to  10  microns 
to  be  representative  of  the  class  of  injectors 
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contemplated  for  this  type  of  device.  Chamber 
pressures  histories  are  shown  for  injection 
durations  of  0,  1,  2  and  3  milliseconds.  As  the 
time  required  for  injection  increases  the  final 
pressure  at  the  end  of  injection  increases. 

Satellite  Applications 

Previous  studies  of  the  potential 
benefits  for  pulsed  combustors,  including  pulse 
detonation  rocket  engines,  have  focused  on 
improvements  in  specific  impulse,  however,  for 
space  applications,  the  improvements  in 
theoretical  performance  are  minimal  because 
large  expansion  ratios  can  be  achieved  with  low 
chamber  pressures  produced  by  pressure-fed  feed 
systems.  Space  payoffs  for  pulsed  rockets  will 
ultimately  be  determined  by  comparing  an 
optimized  pulsed  rocket  system  with  an 
optimized  conventional  system,  however,  in 
order  to  determine  if  there  is  enough  potential  to 
warrant  further  investment,  basic  sensitivity 
studies  such  as  the  following  must  first  be 
performed. 

We  consider  a  typical  application  for 
chemical  propulsion  on  a  satellite:  a  pressure- 
fed,  450  Newton,  apogee  engine.  We  consider 
two  scenarios  for  replacing  the  existing  thruster 
with  a  pulsed  thruster.  In  the  first  scenario  the 


Constant 

Pulsed 

Pulsed 

Pressure 

Engine  - 

Engine 

Baseline 

Scenario 

Scenario 

1 

2 

Thrust  (N) 

450 

450 

450 

Min.  Pc  (MPa) 

0.69 

0.69 

0.26 

Max.  Pc  (MPa) 

0.69 

3.00 

1.10 

Flow  (Kg/s) 

0.143 

0.141 

0.143 

Area  Ratio 

164 

375 

164 

Throat  Dia.  (cm) 

1.91 

1.26 

1.91 

Exit  Dia.  (cm) 

24.5  , 

24.5 

24.5 

Chamb.  Dia. 

(cm) 

5.08 

3.30 

5.08 

Length  (cm) 

60.5 

57.2 

60.5 

Specific  Impulse 

316 

322 

316 

Engine  Mass 
(Kg) 

4.0 

5.3 

4.75 

Tank  Mass  (Kg) 

22.2. 

22.2 

15.9 

Table  4  Notional  study  of  pulse  combustor 
replacements  for  an  apogee  engine 


spacecraft  manufacturer  wishes  to  increase  the 
Isp  of  the  thruster,  but  cannot  accept  an  increase 
in  the  nozzle  exit  area  due  to  space  envelope 
restrictions.  The  manufacturer  wishes  to  retain 
the  existing  propellant  systems  as  much  as 
possible,  including  propellants,  operating 
pressures  and  flowrates.  In  the  second  scenario, 
the  manufacturer  would  like  to  reduce  feed 
pressures  to  reduce  feed  system  component 
weights  without  sacrificing  performance. 

Engine  performance  is  calculated  using 
the  model  described  above.  Dimensions  and 
weights  are  obtained  using  the  PRO:AIAA 
program.  Results  for  the  baseline,  constant- 
pressure  case  and  the  pulsed  engine  scenarios  are 
given  in  the  table. 

Consistent  with  the  stated  requirements, 
the  thrust  is  the  same  in  each  case.  The  space 
envelope  is  nearly  the  same  with  the  nozzle  exit 
diameter  the  same  for  each  engine  and  the 
overall  length  reduced  by  3  cm  for  pulsed  engine 
scenario  1. 

The  goals  for  scenario  1  are  achieved 
through  an  increase  in  the  average  chamber 
pressure,  allowing  the  same  thrust  to  be  achieved 
with  a  smaller  throat,  which  in  turn  allows  a 
higher  expansion  ratio  nozzle  to  be  used  without 
changing  the  exit  diameter.  This  produces  an  Isp 
improvement  of  6  seconds.  The  minimum 
chamber  pressure  was  assumed  to  be  determined 
by  the  feed  system  and  was  set  at  the  same  value 
used  for  the  constant  pressure  engine.  The 
maximum  chamber  pressure  was  determined  by 
a  practical  upper  limit  on  expansion  ratio  of 
approximately  375.  Engine  mass  is  increased 
by  1.3  Kg  as  a  result  of  the  higher  peak  chamber 
pressure. 

The  goals  for  scenario  2  are  achieved  by 
decreasing  the  pressure  of  the  feed  system 
allowing  a  reduction  in  tank  mass  of  6.3  Kg. 

The  average  chamber  pressure  is  the  same  as  the 
baseline  case  producing  an  engine  of  identical 
dimensions  and  performance.  The  increase  in 
peak  pressure  results  in  an  increase  in  engine 
mass  of  0.75  Kg,  resulting  in  an  overall  decrease 
in  mass  of  5.5  Kg. 

The  significance  of  these  changes  must 
be  considered  within  the  context  of  satellite 
economics.  The  useful  lifetime  of  a  satellite  is 
generally  limited  by  the  propulsion  system.  A 
rule  of  thumb  for  satellite  operations  is  that  a  1 
second  improvement  in  Isp  corresponds  to 
approximately  50  days  of  station-keeping  and  a 
six  second  improvement  corresponds  to  an 
additional  year.  Also,  one  month  on  station  in 
geosynchronous  orbit  requires  approximately  2 
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Kg  of  propellant  and  5.5  Kg  of  additional 
propellant  corresponds  to  approximately  three 
additional  months.  Each  additional  month  on 
station  is  worth  several  million  dollars  in 
revenue.  The  potential  payoff  from  the  two 
scenarios  are  in  the  range  of  tens  of  millions  of 
dollars  per  satellite,  fully  justifying  the  cost  of 
developing  these  systems. 


Conclusion 

It  has  been  shown  that  the  potential 
benefits  of  pulse  detonation  propulsion  can  also 
be  realized  using  a  simpler  explosion-cycle 
system.  The  ignition  system  and  buffer  gas 
injection  can  be  eliminated.  Extremely  high 
pressures  associated  with  detonations  of 
condensed  phases  can  be  avoided. 

The  specific  impulse  of  a  pulsating 
rocket  is  the  same  as  a  steady  state  rocket 
operating  with  the  same  average  mass  flow  and 
throat  area.  This  sets  the  value  of  the  equivalent 
steady  state  chamber  pressure  which  determines 
the  performance.  The  nozzle  area  ratios 
optimize  at  the  same  value. 

There  are  three  characteristic  times  that 
control  explosion-cycle  engine  operation: 
injector  pulsing  period,  chamber  blowdown  time, 
and  heat  release  time.  The  heat  release  time 
depends  on  both  the  intrinsic  burning  rate  of  the 
propellant  and  the  time  required  for  injection. 

For  the  cycle  to  operate  the  heat  release  time 
must  be  less  than  the  chamber  blowdown  time, 
which  in  turn  must  be  less  than  the  injector 
pulsing  period.  The  magnitude  of  chamber 
pressure  oscillation  was  shown  to  depend  on  two 
non-dimensional  time  scales,  the  ratio  of 
chamber  blow  down  time  to  the  injector  pulsing 
period,  and  the  ratio  of  the  heat  release  time  to 

the  chamber  blow  down  time. 
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-Figure  L  Effect  of  injection  frequency  on  chamber 
pressure  for  a  fixed  engine  geometry.  Heat  release 
occurs  instantaneously. 


Figure  2.  Specific  impulse  as  a  function  of  ambient 
pressure.  Performance  advantage  of  pulsed  propulsion 
increases  at  high  ambient  pressure. 
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Figure  3.  Chamber  pressure  oscillations  for  chambers 
with  varying  geometry  and  operating  frequency.  The 
magnitude  of  pressure  oscillations  is  a  function  of  the 
blowdown  time  and  injector  pulsing  period  only. 


Figure  4.  Chamber  pressure  plotted  in  non-dimensional 
coordinates  showing  that  the  pressure  history  is  a 
function  of  the  ratio  of  blowdown  time  to  injector 
nul  si  ns  neriod. 
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Figure  5.  Effect  of  varying  the  ratio  of  blowdown 
time  to  injector  pulsing  period  in  the  limit  of  short 
heat  release  time. 


Figure  6.  Effect  of  finite  heat  release  time  on  chamber 
pressure  oscillations 
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Figure  7.  Effect  of  finite  injection  duration  on  rate 
of  chamber  pressure  rise.  Chamber  pressure  must 
not  exceed  injector  discharge  pressure  during  the 
injection. 
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